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A neutral particle with general spin and magnetic moment moving in an 
fT)- arbitrarily varying magnetic field is studied. The time evolution operator for 

the Schrodinger equation can be obtained if one can find a unit vector that 
. satisfies the equation obeyed by the mean of the spin operator. There exist at 

least 2s + 1 cyclic solutions in any time interval. Some particular time interval 
may exist in which all solutions are cyclic. The nonadiabatic geometric phase 
iy) | for cyclic solutions generally contains extra terms in addition to the familiar 

one that is proportional to the solid angle subtended by the closed trace of the 
spin vector. 
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Since the discovery of the geometric phase [3, Q, H, 0, H, S S B 0| 5 particles with spin and 
ma gnet ic moment movin g in time-dependent mag netic fields have received much attention 
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2ll|22L l23j]. though the subject is rather old and some 



discussions can be found in the textbook |24| . Neutral particles are of special interest since the 
problem is easier and the Schrodinger equation can be solved analytically in some special cases, 
say, uniform magnetic fields with a fixed direction or rotating ones. Thus the model is very 
qh| suitable for the study of time evolution, cyclic solutions and geometric phases etc. However, 
some problems in this model are still not clear. First, the Schrodinger equation for general 
k>( ' spin, or even for spin 1/2, in an arbitrarily varying magnetic field seems impossible to be solved 
analytically. Second, though the existence of cyclic solutions in a given time interval may be 
ensured by the existence of eigenvectors for the unitary time evolution operator, it seems not 
clear there are how many ones in the general case. Third, for spin 1/2 it is well known that any 
cyclic solution in an arbitrary magnetic field has a nonadiabatic geometric phase proportional 
to the solid angle subtended by the closed trace of the spin vector. For hig her spin, however, 
this is true only for cyclic solutions with special initial conditions 

HE EES- 

For more general 

cyclic solutions in a rotating magnetic field, we have shown that the nonadiabatic geometric 
phase contains an extra term in addition to the one proportional to the solid angle. The extra 
term vanishes automatically for spin 1/2. For higher spin, however, it depends on the initial 
condition j^. It is still not clear what is the relation between the nonadiabatic geometric 
phase and the solid angle for general cyclic solutions in an arbitrary magnetic field. In this 
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paper we are going to deal with these problems, and try to solve them to some extent. Besides 
the theoretical interest in itself and other applications , this subject has been recently 

recognized to be of great interest in the physics of quantum computation [3, . 

Consider a neutral particle with spin s (s = 1/2,1,3/2,...) and magnetic moment /x = 
fis/s, where s is the spin operator in the unit of h, satisfying [s i5 Sj] = ie^kSk (for spin 1/2, 
s = a/2 and /x = ficr). In a uniform but time-dependent magnetic field B(£) = B(t)n(t) 
where n(£) is a unit vector, it has the Hamiltonian H(t) = [i ■ B(£) = —hu>B{t)s ■ n(t), where 
unif) = fJ>B(t)/sh, and the Schrodinger equation ihd t ip{t) = H{t)ip{t) takes the form 

d t ip(t) =iuj B (t)s-n(t)xp(t). (1) 

Define the spin vector as 

v(t) = (V>(*W(*))- (2) 
Using Eq. (JTJ) it is easy to show that it obeys the equation 

v(t) = -w B (t)n(t) x v(t). (3) 

We are not going to solve Eq. (0) in the general case since this seems impossible. However, 
we will show that the time evolution operator for Eq. (0) can be obtained without any 
chronological product if one can find one nontrivial (nonzero) solution, say, a unit vector eft), 
to Eq. (|3*|) . This is of interest since the latter is easier and more cases can be solved |l4j . 
Actually, Eq. involves operators while Eq. Q involves only c-numbers. On the other 
hand, if transformed to a matrix equation, Eq. (Q) involves 2s + 1 complex variables, while 
Eq. (jHI) involves only three real ones (actually two since it is easy to see that v 2 (t) = v 2 (0)). 
Using the time evolution operator and Eq. (jBJ), one can discuss cyclic solutions and geometric 
phases in a most general way. In particular, we will show that there exist at least 2s + 1 cyclic 
solutions in any time interval. A general relation between the nonadiabatic geometric phase 
and the solid angle subtended by the closed trace of the spin vector will be established. 

To begin, we take an arbitrary unit vector e , and the eigenstate of s • e with eigenvalue 
m will be denoted by Xm- We take the initial state of the system to be 0(0) = Xm, that is 

s • eoV'(O) = 777/0(0), m — s, s — 1, . . . , — s. (4) 

Obviously v(0) = meo in this initial state. Now we define a vector e(t) by Eq. (jH)) with the 
initial condition eo, that is 

e(t) = -co B (t)n(t) x e(t) (5) 

with e(0) = e . We would assume that B(t) varies continuously, so that any solution e(t) is 
well behaved. As pointed out above, e 2 (t) = e^, so e(t) is a unit vector at any time. We have 
proven in Ref. [22 1 that 

s ■ e(t)if>(t) = mip(t) (6) 

holds at all later times. To be self-contained, we repeat here the proof by induction. 

By definition, Eq. (JHJ) is valid at t = 0. We assume that it is valid at time t, what we need 
to do is to show that it is also true at time t + At where At is an infinitesimal increment of 
time. In fact, using Eqs. ((TJ) and (0) we have 

i/j(t + At) = ip(t) + iuj B (t)s ■ n(t)ip(t)At, (7a) 



e(t + At) = e(t) - u B (t)n(t) x e(t)At. 



(7b) 



After some simple algebra, the conclusion is achieved. 

Because e(t) is a unit vector, we can write in some rectangular coordinates 

e(t) = (sin#(t)cos0(t), sin 9(t) sin <j>(t), cos9(t)). (8) 

Using the formula 0, 

e i 5s .b se -i ?s .b = [ s _ ( s . b)b] cos £ + (b x s) sin f + (s • b)b, (9) 

where b is any unit vector, it is not difficult to show that 

s-e(t) =e -ms<t) Sze w(t)s- d (t)^ (10) 

where 

d(t) = (-sin0(t),cos0(t),O). (11) 
Therefore the eigenstate of s ■ e(t) with eigenvalue m is 

where Xm is the eigenstate of s z with eigenvalue m, and a m (t) is a phase that cannot be deter- 
mined by the eigenvalue equation. However, a m (t) is not arbitrary. To satisfy the Schrodinger 
equation, it should be determined by the other variables 6{t) and 4>{t). In fact, the above 
equation yields 

(V(t), V(* + At)) = 1 + id m (t)At + e ie W s - d Wa t e- ie W s - d W Xr °J At. (13) 
Using the formula [0 

e- p W^e Ji 'W= C e- XF ®F(t)e XF ® dX, (14) 

where F(t) is any operator depending on t, and then using Eq. (JHJ), we obtain 

(ip(t),ip(t + At)) = l + id m (t)At + im[l -cos0(t)]0(t)At. (15) 

On the other hand, from Eq. (0) we have 

(V>(t), ^(t + At)) = 1 + iw B (t)v(t) • n(t)At. (16) 

Note that v(t) and e(t) satisfy the same equation, and v(0) = me , we have v(t) = me(t). 
Comparing the two results above and taking this relation into account, we obtain 

a m (t) = -m[l - cos9(t)]i>(t) + muj B {t)e{t) ■ n(t). (17) 

Therefore 

a m (t) - a m (0) = ma(t), (18) 

where t 

a(t) = - [ [1 - cos0(t')]0(t') dt' + f u B {t')e(t') ■ n(t') dt'. (19) 
Jo Jo 

Substituting into Eq. fH2|) we obtain 

^(t) = e-w(*)-d(*)eia(t)-. e itf(o)-d(o)^( 0) _ (20) 



We denote the time evolution operator as U(t), denned by the equation ip(t) = U(t)ij)(Q) with 
an arbitrary "0(0), then the above equation is equivalent to 



U(t)x 



m 



e -ifl(i)B.d(t) e ia(t)a, e W(0)B.d(0) 



Xm- 



(21) 



Now an arbitrary initial state ip(0) can always be expanded as 



V>(0) = ^CmXm- 



(22) 



ni 



Applying U(t) to both sides of this equation, using Eq. (|2*T|). and noting that the operators 
on the right-hand side of that equation is independent of m, we immediately realize that Eq. 
(}20j) is in fact valid for an arbitrary initial state. Thus we arrive at the result 



Eq. (|23bj) is suitable for the general discussions below while Eq. ()23aj) may be more convenient 
for practical calculations. 

Let us make some remarks on the result. First, we see that once a solution of Eq. (jSJ) is 
found, the time evolution operator for Eq. (P) is available and it involves no chronological 
product. The result depends formally on eo, but eo is merely an auxiliary object, hence the 
result must be essentially independent of it, though it might be difficult to prove this explicitly. 
In practical calculations, one should choose a solution e(t) that is as simple as possible such 
that U(t) can be easily reduced to the simplest form. When this approach is used to the 
simple cases such as rotating magnetic fields or ones with a fixed direction, it indeed leads to 
the same results as those obtained previously. Second, the operator U(t) depends not only 
on e(t), but also on the history of it. This is obvious from Eq. |T9|) . Third, though (f)(t) is 
indefinite when 6(t) = or tt, the above result U(t) is in fact well behaved everywhere. That 
it is well defined at 8(t) = is obvious. If 6(0) = 0, there is no problem either. The case with 
#(0) = tt can be avoided since one can always choose a coordinate system such that 8(0) ^ tt. 
However, for a general evolution, the case with 9(t) = tt at some instant cannot be avoided. 
Thus we must show that U(t) is well behaved at 9(t) = tt. Suppose that 6(to) = tt, then we 
have 0(td) — 0(^o ) = ±7r, d(tg ) = — d(tQ ), and a(t~Q ) — a(t^ ) = ^2tt. With these relations 
it is not difficult to show that U(t$) = U(t ). Since both U(t^ ) and U(t$) are well defined, 
we may define U(to) = lim^ <0 U(t). This makes U(t) well defined and continuous at t = to. 
Fourth, by straightforward calculations it can be shown that dtU(t) = iu>B(t)s ■ n(t)U(t) and 
U(0) = 1, as expected. 

Now we can go further to discuss cyclic solutions in any time interval [0, r] where r is an 
arbitrarily given time. These cyclic solutions are not necessarily cyclic in subsequent time 
intervals with the same length, say, [r, 2r] . 

Since Eq. (jSJ) is a linear differential equation, the general solution e(t) must depend on the 
initial vector eo linearly. Thus it can be written in a matrix form 



U(t) = e 



i0(t)s-d(t) e ia(t)s z e i6»(0)s-d(0) 



(23a) 



Using Eq. Q, it can be recast in the form 



U(t) = e 




(23b) 



ei(t) = Eij^eoj, 



(24) 



where the matrix E(t) is obviously real. If both e\(t) and ^(t) are solutions to Eq. (0), it is 
easy to show that e\(t) ■ e2(t) = ei(0) • e2(0). Therefore the matrix E(t) is an orthogonal one, 



and its eigenvalues at any time t has the form {1, cr(t), a*(t)}, where a(t) is a complex number 
with \a(t) \ = 1, and cr*(t) its complex conjugate. 

If cr(r) 7^ 1, one eigenvector 77 (r) of the matrix -E(r) with eigenvalue 1 can be found, which 
satisfies £'y(r)^-(r) = ^(t). It can be taken as real and normalized. Now if we choose 

e = 77(r), (25) 

we have e^r) = Eij(r)e j = Eij(T)r}j(r) = ^(t) = ej(0), that is 

e(r) = e . (26) 

This means that 9(t) = 6(0) and d(r) = d(0), and leads to 

C/( T ) = e iQ ( T ) s - eo . (27) 

Now it is clear that with the initial condition ijj(0) = Xm (m = s,s — 1, . . . , — s), we have a 
cyclic solution in the time interval [0,r]. More specifically, ipij) = e" 5 V ; (0), where the total 
phase change is 5 = ma(r), mod 2ir, with a(r) given by 

a(r) = -fi e + / w B (*)e(t) • n(t) dt, (28) 



where 



fi e = / [l-cos0(t)]0(t)dt 



is the solid angle subtended by the closed trace of e(t). Notice that v(£) = me(t), the dynamic 
phase (3 = —hr 1 JJ" (H(t)) dt turns out to be 

P = m f uj B (t)e(t) ■ n(t) dt. (29) 
Jo 

Therefore the nonadiabatic geometric phase is 

7 = 5-/3 = — mQ e , mod 2tt. (30) 

Since Q e = e(m)n„,mod 47r, where Q v is the solid angle subtended by the closed trace of the 
spin vector, we have finally 

7 = — |m|f2„, mod 2tc, (31) 

in accord with the results previously obtained Thus we see that 2s + 1 cyclic 

solutions are available in any time interval [0,r], and all phases can be expressed in terms of 
the unit vector e(i). 

States with initial condition other than the above ones are in general not cyclic ones, even 
those in which v(0) points in the direction of e such that v(r) = v(0). However, if a{r)/n 
happens to be a rational number other than an even integer (the case with a(r)/7r an even 
integer will be discussed below), some other cyclic solutions may be available. To be more 
specific, let a(r) = pn/n, where n is a natural number and p an integer. When n = 1, p is an 
odd number, and when n > 1 it is prime with p. If s > n, we have cyclic solutions with initial 
condition, say (no such solution exits for s = 1/2), 



N 2 

^(0) = C i*m+2ni, 



(32) 



A ? 2 



where Ni and N 2 are nonnegative integers, and 

m - 2nN 1 > -s, m + 2nN 2 < s, V] |c,-| 2 = 1. 

j = -ATl 

In this initial state 

v(0) = v e , (33) 

where v = ^2f=_ Nl (^n + 2nj)\cj\ 2 may be either positive or negative. It is easy to show that 
V>( T ) = e "V(0)> where 5 = ma(r),mod 27r. Because of Eq. (}33|) . we have v(t) = t>oe(t), and 
the dynamic phase is 

PT PT PT 

/3 = -h~ 1 (H(t))dt = uj B (t)v(t) ■ n(t) dt = v / w fl (t)e(*) ■ n(t) dt. 
Jo Jo Jo 

Using Eq. (1281). we have 



/3 = i;o[a(r) + n e ]. (34) 

This holds regardless of the values of vo and a(r), as long as Eq. ((33|) is valid. Suppose that in 
the process from t = to t = r, e(t) encircles the polar axis K times (K > for anticlockwise 
traces and K < for clockwise ones), then we have 

Q e + fi_ e = 4nK. (35) 

This leads to fo^e = |fo|^u + 2ttK(vq — \vq\). The geometric phase turns out to be 

7 = — l^ol^ + 27ri^(|f | — v ) + (m — Vo)pii/n, mod 2n. (36) 

In this case we see that the geometric phase contains extra terms in addition to the one 
proportional to Q v , unless the sum of these extra terms happens to be an integral multiple 
of 2tt. Note that the above relation holds for the case v o = as well, though Q v is not well 
defined in this case. 

If it happens that a(r) = 2krc, then U(t) = e l2nks becomes a c-number and all solutions 
are cyclic in the time interval [0, r]. However, this is true only on the premise of (|26jl . Thus 
a(r) = 2kn alone is not a sufficient condition for all solutions to be cyclic, but it is easy to 
show that it is a necessary one. A sufficient condition is <t(t) = 1. 

Now if er(r) = 1, E{r) becomes a unit matrix. In this case any vector is its eigenvector 
with eigenvalue 1. Therefore Eqs. (f2T)|) and (|27j) hold for any unit vector eo- In particular, we 
have 

JJ(r) = e iai ^ Sx = Q ia2( ~ T ^ s y = e iQ 3(r)s z ^ 

where a>i{r) is given by Eq. (|28|) with e = e x and similarly for 02(1") and 013(1"). Since s x , s y 
and s z are independent operators, the above equation cannot hold unless a«(r) = 2nki where 
ki are integers such that U{r) becomes a c-number. In general, with an initial unit vector eo, 
we have 

a{r) = 2nk, (38) 

and 

U(t) = e ia(T)s - eo = e i2 " fcs . (39) 

Let us see what is the dependence of a(r) or k on the direction of e . Consider two initial 
unit vectors e and e , whose difference Se = e' Q — e is infinitesimal (then 5e • e = 0). The 
difference in a(r) is, according to Eq. ((21 



Sa(r) = o!{t) - a(r) = -5il e + / u B (t)6e(t) ■ n(t) dt. (40) 

Jo 



Since e^t) = Eij(t)e j, we have 5ei(t) = Eij(t)8e j, and the second term in the above equation 
is consequently infinitesimal. Moreover, the trace of e'(i) is very close to that of e(i), thus the 
difference in the solid angles subtended by them is infinitesimal. Therefore both 5a(r) and 
5k = 5a(r)/27i are infinitesimal as well. Now that k can take only integer values, an obvious 
consequence is that 5a{r) = and k! = k. 

A subtle case has been overlooked in the above discussions, however. This happens when 
e(t) goes by the south pole {6 = n) on one side and e'(t) goes by it on the other. In other 
words, the closed trace of e(t) as a whole goes across the south pole when e varies to eg. In 
this case a(r) will be changed by an integral multiple of An and k by an even integer. One 
can of course removing this change by rotating the coordinate system such that the trace of 
e(t) does not go across the south pole. However, this can only be done locally, and globally it 
is impossible in general. In other words, one cannot choose a coordinate system such that k is 
the same for all e , except for some simple cases, say, magnetic fields with a fixed direction. 
The above dependence of «(t) on eo has no consequence on the time evolution operator as 
expected. This is easily seen from Eq. when k changes by an even integer, U(r) remains 

the same. 

In the special case where k is the same for all eo, a geometric explanation of k is available. 
If we take the unit vector e' = — eo as an initial condition to Eq. (|5|) . then the solution is 
e'(t) = — e(t). Since a(r) corresponding to e(t) and a'(r) corresponding to e'(t) are equal as 
assumed, we have, according to Eq. (J2EJ), 2a (r) = — (f2 e + f2_ e ). As before, Q e + f2_ e = 4iiK, 
and 2a{r) = 47r/c, so that k = —K. Here K is the winding number of e(t) around the polar 
axis. Unfortunately, the geometric meaning of k in the general case is still not clear. 

Now that U(t) is a c-number, all solutions to Eq. become cyclic in the time interval 
[0,r]. The total phase change is 5 = 2nks = sa(r),mod 2ir. Since this phase is determined 
only up to an integral multiple of 2ir, the dependence of a(r) or k on the initial vector eo does 
not affect the result. For convenience we take the a(r) or k that is the same as that appears 
below. If in the initial state v(0) ^ 0, we take e = v(0)/t> where v = |v(0)| > 0. Thus 
v(t) = foe(t) and fl e = fl v . The dynamic phase is, as shown in Eq. (JMj) . (3 = v [a(r) + fl v ]. 
It should be remarked that both a(r) and Q v may be different by an integral multiple of 4-7T 
in different coordinate systems. However, (3 has a definite value, independent of the choice of 
coordinate systems, which can be easily seen from the expression (3 = v JJ" u>B(t)e(t) ■ n(t) dt. 
The geometric phase turns out to be 

7 = — v Q v + (s — Vo)2iik, mod 2n. (41) 

Here the first term is the familiar one, but an extra term appears. If s = 1/2, we have Vq = 1/2 
for any initial state, and the above result reduces to 7 = —Q v /2, a well-known result. For 
higher spin, however, the extra term depends on the initial state, as both vq and k depend on 
it. It vanishes (mod 2n of course) when s — t>o is an integer, especially when the initial state 
is an eigenstate of s • e (it cannot be an eigenstate of s ■ e with some other unit vector e' 
since v(0) points in the direction of e ), as expected. If in the initial state v(0) = 0, we have 
v(t) = and thus (3 = 0. It is easy to see that Eq. (|4Tj) holds in this case as well, though Q v 
is not well defined. 

The general result (J4~T|) has been confirmed by practical calculations in the case of a rotating 
magnetic field, where both 7 and Q v can be calculated explicitly j22j. For a magnetic field 
with a fixed direction, it is easier to carry out similar calculations to verify this result. 

In summary, we have shown that the time evolution operator for the Schrodinger equation 
(JTJ) can be obtained if one nontrivial solution to Eq. (jSJ) can be found. We proved that at 
least 2s + 1 cyclic solutions of the Schrodinger equation exist in any time interval. These cyclic 



solutions can be worked out in principle if the general solution to Eq. (J3J) is known. There 
may exist some particular time interval in which all solutions are cyclic. The nonadiabatic 
geometric phase for cyclic solutions contains in general extra terms in addition to the familiar 
one that is proportional to the solid angle subtended by the trace of the spin vector. For spin 
1/2 there is no such extra term. 
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